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I Abstract. In this paper we discuss the consequences of a Killing symmetry on the 

■ local geometrical structure of four-dimensional spacetimes. We have adopted the point 
. of view introduced in recent works where the exterior derivative of the Killing plays a 

00 ' fundamental role. Then, we study some issues related with this approach and clarify 

why in many circumstances its use has advantages with respect to other approaches. 
, We also extend the formalism developed in the case of vacuum spacetimes to the general 

' case of an arbitrary energy-momentum content. Finally, we illustrate our framework 

C^j . with the case of spacetimes with a gravitating electromagnetic field. 
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1. Introduction 

The imposition of symmetries is and has been one the most successful ways of simplifying 
the theory in order to deal with the essence of physical problems. In the case of General 
Relativity they have been extensively used to simplify Einstein's field equations, and 
in this way to obtain exact solutions describing the gravitational field of the idealized 
situation. Of course, due to the non-linear character of the field equations, assuming 
many symmetries can lead to solutions which may not capture the actual behaviour of 
the physical system they are meant to represent. Then, it is possible that deviations 
from the idealized situation behave in a very different way. Nevertheless, we know of 
many systems whose dynamics is driven by the gravitational interaction and where there 
is a symmetry that is present in the spacetime region of interest. Then, it is important 
to know how the imposition of a symmetry affects to the structure of the spacetime and 
which can be the consequences for a specific problem. 

In this paper we will consider spacetimes with only one Killing symmetry, which 
can describe many physical situations of interest. Examples of present interest are 
provided by physical systems in quasi-equilibrium configurations, like neutron-star 
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binaries around the innermost stable circular orbit (ISCO), which are of great interest 
in the study of gravitational wave emission by binary systems. In this context, methods 
to construct such configurations in order to locate the ISCO rely on the assumption 
of the existence of an approximate Killing vector M. Other problems of interest in 
which a symmetry is present are those involving axisymmetric configurations, stationary 
systems, etc. 

There are many papers in the literature dealing with spacetimes possesing Killing 
vector fields (KVFs hereafter). Most of these papers are devoted to the search of 
exact solutions of Einstein's field equations, and KVFs were considered as a way of 
simplifying the problem (see f2j). Then, the main aim was to develop and apply 
techniques to solve the Einstein equations when a group of Killing symmetries was 
imposed. Despite the huge amount of work done there are few studies that concentrate 
only on the study of spacetimes containing just one KVF. Then, few is known about 
the geometrical properties of a Killing symmetry in connection with the spacetime 
geometry. Examples of works where formalisms for one KVF are introduced are the 
works by Collinson and French || (for conformal KVFs) and by Perjes [|J. In the first 
work the equations (including the Killing equations) are written using the Newman- 
Penrose formalism ||, which leads to simplifications when the NP basis is adapted to 
the algebraic structure of the spacetime. However, it has the limitation that it does not 
introduce any quantity adapted to the (conformal) KVF. Then, by one hand it can be 
used for groups of (conformal) KVFs but, on the other hand, it cannot deal with the 
particular characteristics of the symmetry. In contrast, the work by Perjes combines 
in a powerful way the spinor formalism (the NP formalism is a byproduct of it) and 
quantities describing the structure of the symmetry. However, since it uses quantities 
related to the Ernst potential, the generalization to general spacetimes with one KVF 
is not obvious at all. 

Recently, we introduced || |7j a new approach to the study of vacuum spacetimes 
with one KVF which is adapted to the structure of the Killing symmetry and which can 
be easily generalized to any matter content (this is one of the subjects of the present 
paper). As we showed in the previous works 0, [7|], this approach naturally leads to a 
classification |7|] of the spacetimes with at least one symmetry, which is based on the 
properties of the KVF and on the relations of its algebraic structure with that of the 
spacetime (the algebraic structure of the Weyl tensor). Moreover, in 0, |7| we showed 
the potential of this approach by proving a number of new results in the vacuum case 
(see ]n| for an error in one of them). Generalizations of this formalism have recently 
appeared for the case of homotheties [pl| and conformal KVFs [[16 . 



The approach introduced in |], 0] considers the exterior derivative of the Killing, 
which we will call the Papapetrou field after its introduction in || (sometimes also 
called the Killing 2- form or the Killing bivector |9|, |Kj), and its algebraic structure, 
as main objects of the study. The Papapetrou field has been used before in the 
search of exact solutions (see, e.g., [|TT]]) and in the study of external magnetic fields 
in black holes [T^]. As we discussed in detail in f| |7|, one main advantage of this 
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approach is that it provides a framework in which one can establish, in a natural way, 
connections between the particular characteristics of the symmetry (the structure of 
the Papapetrou field) and the algebraic structure of the spacetime (the Petrov type 
of the Weyl tensor which, in particular, leads to a detailed classifications of the 



spacetimes with an isometry 0. Moreover, with the introduction of the Papapetrou field 
as a main quantity we constructed a new formalism for the study of vacuum spacetimes 
with an KVF [0. It is an extension of the Newman-Penrose formalism ||, where new 
variables associated with the KVF and its Papapetrou field are introduced. Within this 
formalisms, the integrability conditions for the components of the KVF are analyzed 
from a new perspective. Using this new point of view it was shown in @ that from the 
integrability conditions one gets expressions for all the Weyl tensor components in terms 
of the other variables, reducing in this way the number of equations to be considered 
(we avoid the use of the second Bianchi identities). 

In this paper we analyze in detail the structure of this approach and the reasons 
why it has advantages over other approaches, and we extend the formalism to spacetimes 
with a KVF and with a generic energy- momentum content. The plan of the paper is 
the following: In section |2| we present the main ideas and justify why the adoption 
of the Papapetrou field as a main quantity means to have advantages when dealing 
with spacetimes with a KVF. At the same time we study in detail, from a new 
perspective, the integrability conditions for the components of the KVF. This part will 
also involve an explanation of previous results 0, in a broader context. In section [3] 
we study the integrability conditions for the Papapetrou field, which satisfies Maxwell 
equations. In section ^ we extend the framework introduced in [0] for vacuum to general 
spacetimes. For the sake of clarity we will give the equations that come out from this 
formalism in the |Appendix A| , where they have been expressed using the Newman- 



Penrose formalism ||. In section |5] we illustrate these developments with the case of 
spacetimes with a gravitating electromagnetic field, and in particular, with the case of 
the Kerr-Newman metric. We finish with some remarks and conclusions in section 
Throughout this paper we adopt units in which c = 8nG = 1 , and we will use the 
conventions and definitions introduced in unless stated otherwise. 

2. On the consequences of the integrability conditions for the Killing 
equations 

From now on we will assume the spacetime is endowed with a non-null KVF £ 

N= gab ce^o. (i) 

From the definition of a Killing symmetry it follows that £ satisfies the Killing equations 

£ ( g ab = ^ £«;6 + & ;a = 0, (2) 

where the semicolon denotes covariant differentiation. We can know about the Killing 
symmetries that a given spacetime admits by solving these equations. Of course, there 
are spacetimes that do not admit any Killing symmetry, in which case one finds that 
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the Killing equations have no solutions. We can investigate whether or not the Killing 
equations admit solutions in a given spacetime by studying their integrability conditions, 
which can be expressed in the following compact form 

£^ be = £$R bed = £ (,R bcd;ei = £$R bcd\e\e2 — " " " — , (3) 

where <£ € denotes Lie differentiation along £, T a bc are the Christoffel symbols, and R a bcd 
the components of the Riemann curvature tensor. As we can see, the integrability 
conditions involve, in general, derivatives of the curvature tensor. 

In recent works || |7| we introduced an alternative way of dealing with vacuum 
spacetimes admitting a Killing symmetry. The main object in this approach is the 
exterior derivative of the KVF 

F = d$, F a b = £,b;a ~ £,a;b = 2£fe ;a , 

where we have used the Killing equations @. This quantity is important for several 

reasons. As it was firstly recognized by Papapetrou [§], a KVF £ can always be seen 

as the vector potential generating an electromagnetic field, F a b, satisfying Maxwell 
equations 

F [ab . c] =0, F ab , b = J a , (4) 

where J a is the conserved current given by 

J a = 2R\e => J a . a = 0, (5) 

being R a b = R c a cb the Ricci tensor. Apart from the Maxwell equations, F a b satisfies the 
following relation 

£ € F<* = 0, (6) 

which is a direct consequence of the Killing equations and it is equivalent to 

eF ab ;e = 0. (7) 

In the approach proposed in || [7[ the components of the Papapetrou field F ab are 
promoted to the level of variables of the formalism, like the components of the KVF. 
Then, taking into account the Killing equations (fj), the equations for £ are now 

&>;a = \F a b with F ab = F[ ab ] . (8) 

In the Appendix A| we have expressed them using the NP formalism. 

Their integrability conditions are equivalent to the Ricci identities for the KVF £ 

£,a;dc £,a;cd Rabcd^, • (9) 

In the four- dimensional spacetimes of general relativity, and taking into account only 
the antisymmetry of the Riemann tensor in the first and second pair of indices 
(Rabcd = R[ab]cd = R a b[cd]) , expression (|) contains 24 independent equations. They 
can be split into four differentiated groups: (i) First Bianchi identities: R{ a bc\d£, d = 
=>• C[ab;c] = 0. They contain four independent equations which lead to the first group 
of Maxwell equations ([|). (ii) The trace: £ a 6 ;c — £ a a] b = Rbd£, d - Four independent 
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components equivalent to the second set of Maxwell equations (f|). (iii) R a bcd£, a £, b = 
=>- C, a (Ca t dc — ia-,cd) = 0. This expression contains 6 independent equations, which 
in combination with those in the group (i), lead to the property (||,[/]) of F a b- (iv) 
The rest, 10 independent equations. Introducing (||) for the derivatives of £ and the 
decomposition of the Riemann tensor into its irreducible parts under the Lorentz group: 
the Weyl tensor C a b c d, the curvature tensor part not determined locally by the energy- 
momentum distribution, the traceless Ricci tensor S a b = R a b — |g a 6-^ an d the scalar 
curvature R = g ab R a b, 

Rabcd = Cabcd + ga[cSd]b - gb[ c Sd]a + J2 (SacSbd ~ SadSbc) R 5 

the remaining 10 equations in (H) can be written in the following form 

C a bcd£, d = —\F a b;c + C[aRb]c + ^ Rd[agb]c ~ |£[ag&]c^ • (10) 

The right-hand side is made out of a term constructed only from the Papapetrou field 
and terms constructed from the Ricci tensor. The last ones can be written in terms 
of the matter variables, described by the energy-momentum tensor T ab , through the 
Einstein field equations 

Gab = Rab — \&abR = Tab , (H) 

where G a b denotes the Einstein tensor. Then, we can think of equation ( |T0"D as a relation 
that will provide expressions for some components of the Weyl tensor in terms of F a b and 
T ab . As we have said before, equation ( TO ) contains 10 independent equations, the same 
number of independent components of the Weyl tensor. Hence, we can determine all 
the components of this tensor from ([3~0|) . However, this is true only in four- dimensional 
spacetimes. For D-dimensional spacetimes we would get an equation analogous to (|I0D, 
with the same left-hand side and with the same terms on the right-hand side but with 
different numerical coefficients (see, e.g., []17| to obtain these coefficients). The number 
of independent components of the Weyl tensor in D dimensions is (see, e.g., 

C(D) = Y2,D(D + 1)(D + 2)(D — 3) (for D > 3) . 

On the other hand, to find the number of independent components in the D-dimensional 
version of equation (pT0|) we can write it as C a b c d^ d = A abc . Then, taking into account 
the information coming from (i)-(iii), the tensor A a b c satisfy the following properties 

A[ab C ] = , A[ ab ]c = A abc , A abc ^ c = , A ab a = . (12) 

From these relations and the fact that we are considering KVFs with a non-vanishing 
norm N [see Eq. (|I])], the number of independent components contained in Eq. (P~0| ) are 

I(D) = iD(D + l)(2D-5). 

Hence, the number of independent components of the Weyl tensor that cannot be 
determined from the integrability conditions ([II]) is given by 

A(D) = C{D) - 1(D) = j^D(D 2 - 1)(D - 4) . 
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Therefore, in the case of General Relativity, D = 4, we can determine completely the 
Weyl tensor [A (4) = 0] in terms of the Papapetrou field F ab and the energy- momentum 
content [see Eq. fllPP]. In higher dimensions!]] this is not true, and we can only determine 
a subset of the independent components of the Weyl tensor. 

As a conclusion we can say that from the integrability conditions flit] ) and in D = 4 
we can find an expression for the Weyl tensor in terms of the Papapetrou field and the 
energy- momentum content. We can find the explicit expression of the Weyl tensor by 
different means. Here we will use a procedure based on the fact that for a given non-null 
vector field we can associate with it what are called the electric and magnetic parts of 
the Weyl tensor. We are interested in the electric and magnetic parts associated with 
the KVF £ : 

^a&[£] — C ac bdC^ d 

Hoife] = *C acbd ^ d (*Cabcd = lVab ef C cdef ) , (13) 

where * denotes the dual operation and r}abcd the components of the completely 
antisymmetric volume four form. The two tensors ( |T3"D are symmetric, trace-free and 
orthogonal to the vector field they are associated with, the KVF in our case, 

£(«&) [t] = Eabfc] , H {ab) [£] = H^t] , E\[t] = H\[i] = , E ab [£]e = H ab [t]Z h = . 

Moreover, they contain all the information about the Weyl tensor. Indeed, each part 
has five independent components and the Weyl tensor has ten. We can reconstruct the 
Weyl tensor from the electric and magnetic parts and the KVF £ through the following 
expression 

C abcd = -^Ct, 9 {(gahef&cdgh ~ VabefVcdgh) E fh [t] ~ (gabefVcdgh + Vabefgcdgh)H fh [t]} , (14) 

where g abcd = g ac g bd - g ad g bc . 

From the integrability conditions ( [T0| ) we can compute the electric and magnetic 
parts associated with the KVF the result is 

E ab [£] = -\¥g ea cd A cdb , (15) 

H ab [t] =-^ e Vea Cd A cdb , (16) 

where the tensor A a b c was introduced above and has the properties shown in ([12]). In 
terms of F ab and T afc it has the following expression 

Aabc = -\Fab;c + Z,[aT b ]c + £ d T d [ a g b ] c ~ f £[agf,] c ^ , (17) 

where T = g ab T ab . Then, introducing this expression in equations (15], Hp and finally 
into the Weyl tensor [Eq. (|T4"D1 we get 

C ab cd = -| { A ab [c ti d] + A cd [a £ b] + 2f A e l c [a A] } • ( 18 ) 

And this is how the integrability conditions for the components of the KVF supply 
an expression for the Weyl tensor in terms of the KVF, the Papapetrou field and the 
energy-momentum tensor. 

|| For D < 3 the Weyl tensor is identically zero. 
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We can insert ([T^) into the contracted second Bianchi identities (which in D = 4 
contain the same information as the non-contracted ones) 



C a bcd' d — Rc[a;b] ~ |g c [a-^,&] ■ 



The result, using Einstein's equations (|iT|) , is 

F a b;c C + Cabcd,F Cd + 2 J[ a;6 ] + \T F ab = . (19) 

One can reach this result starting from Maxwell equations (|J), taking the covariant 
derivative of F[ ab . c } = and using the Ricci identities. This means that fllTf ) is an 
identity as long as F ab satisfies Maxwell equations. 

3. On the Integrability conditions for the Papapetrou field and other issues 

In the previous section we have studied the integrability conditions for the components 
of the KVF. The next step is to study the equations for the Papapetrou field F ab , i.e., 
Maxwell equations (|j). Their integrability conditions have been studied in great detail 
in Here, we will apply their results to the case of a Papapetrou field. To that end, 
it is very convenient to consider the following two points. First, instead of using F ab , we 
will work with the self-dual Papapetrou field 

Fab = F ab + i* F ab (*F ab = \n a bc d F cd ) . 

This will lead to more compact expression in our study. In particular, the Maxwell 
equations look simply as 

F ab - b = J a , (20) 

where J a are the conserved sources, given in @. The second point is to take into 
account the algebraic structure of the Papapetrou field, which in || [jj has been shown 
to be fundamental. The idea is to distinguish between the two possible algebraic types 
of F a b. (i) The regular type, characterized by F ab F a b ^ . In this case we can pick a 
Newman-Penrose basis {k,£,m,fh} so that F takes the following canonical form 

F ab = TW ab , where T = -(^ + i0) and W ab = -2k [a £ b] + 2m [a m b] , (21) 

and where and j8 are the real eigenvalues of F ab , and k and £ are its null eigenvectors 
(principal null directions), (ii) The singular type, characterized by F ab F a b = 0. Now, 
we can choose the Newman-Penrose basis so that F can be cast in the form 

Fab = ZOVab , where V ab = 2k[ a m b ] , (22) 

where 9 is a complex scalar and k gives the only principal null direction. 

The NP basis in ( PT|P2"D are not completely fixed. The transformations that keep 
F a b m the canonical form were discussed in 0j. We only mention the fact that whereas 
in the regular case the eigenvalues ^ and j8 are invariant under these transformation, 
in the singular case the scalar 6 is not, actually we can scale it in an arbitrary way. 

We can now study the integrability conditions of the Maxwell equations for the 
Papapetrou field (p0|) . We are going to consider each algebraic case separately. 
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3.1. Regular case 

Taking into account the following property of W ab 

W a b W bc = g ac , (23) 
one can show that the Maxwell equations (pO|) are equivalent to the following equations 



T,« = Th a + q a , (24) 

where 

K = W ab , c W bc , q a = W ab J b = 2W a b RU c . (25) 

Remarkably, Maxwell equations provide an expression for all the derivatives of the 
eigenvalues $ and (3 (^4|) [their expressions in the NP formalism are given in the 
Appendix A| , Eqs. ([A.11|4A~7[4"|) 1. In this case the integrability conditions for the Maxwell 



equations reduce to the integrability conditions for the complex scalar T . They will 
impose conditions on h a and q a . Following ]nj|, we divide the problem into two cases: 
(a) dh = =>- h^h] = 0. (b) h[ a<b ] ^ 0. Case (a) represents a particular situation in 
which a condition on the 2-form W ab has been imposed. The integrability conditions for 
the complex scalar T are simply 

dq + q A h = q [ci)b] + q [a h b] = . (26) 

This condition is automatically satisfied when q a = 0, or equivalently when R ab ^ b = 
(the KVF is an eigenvector of the Ricci tensor with zero eigenvalue), which includes 
the vacuum case. The integrability conditions (|26f) are first-order equations for q and 
their integrability conditions, d 2 q = , are identically satisfied by virtue of dh = . 



The expressions of this last set of equations in the NP formalism is given in [Appendix 



[A], Eqs. ( |A.34ffA~38| ). Case (b) constitutes the generic situation. The integrability 



conditions for T can be written as follows (see [|19] for details) 

Qabh[ C: d];e ~ ^[a,b] [^cd;e ~ h c Vt de + qjl[d,e]} = , (27) 

where 

Qab = q[a,b] + q\ah b ] ■ 

Another important issue that we must consider is the analysis of the consequences 
of the relation ([?[). To that end it is useful to use the complex 2-forms {U, V, W} , 
where the 2-form U is given by 

U a b = —2l[ a fhb] , 

and W and V have been introduced in (^l|) and (^) respectively. Using (|24]), the 
relation (0) leads to the following three equations 

U ab CW ab . tC = 0, (28) 
V ab CW ab . c = 0, (29) 
C(Th a + q a ) = 0^ CT, a = . (30) 
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These relations mix algebraically the components of the KVF with components of the 
connection and the Ricci tensor. We have given their explicit form in the NP formalism 
in |Appendix A] [Eqs. ( |7n>9T]A~3Tl) 1. 



Finally, from ( Jl5| , p!6| ) and flT6| ) one can construct the invariants of the Weyl tensors 
in terms of the Papapetrou field and the energy-momentum tensor. In the case of 
vacuum spacetimes and a timelike KVFs, one can also construct the Bel-Robinson 
superenergy (see p0| ^T[1 ) associated with the observers following the trajectories of 
the KVF. These quantities have been used in numerical simulations to track, e.g., 
the propagation of gravitational waves. In the regular case and vacuum, by virtue 
of (p4|), the Bel-Robinson superenergy has the form: + {3 2 ) 2 (Connection terms) 2 . 
For the nonvacuum case, the Bel-Robinson tensor is not longer conserved, but there 
are generalizations. In the case of spacetimes with vanishing scalar curvature, which 



include Einstein-Maxwell spacetimes, we have a generalization which also has a 
positive timelike component (superenergy), where apart from terms of the form showed 
before we have terms constructed from the energy-momentum tensor and mixed terms 
of the form: (jz£ 2 + f3 2 ) (Connection terms) (Components of T ab ). 

3.2. Singular case 

The situation in the singular case is completely different. The main reason is that V a b is 
a singular 2- form which obviously does not satisfy (f2"3[). The main consequence of this 
fact is that the Maxwell equations (^0|), which now look as follows 

V a b e fi = -9V a b , b + \J a . (31) 

does not provide expressions for all the derivatives of the complex scalar 9 [see [Appendix] 



A.2|| . Therefore, we can only talk about the compatibility conditions for the expressions 
of the known derivatives of 9 . Actually, as we have previously said, contrary to what 
happens with T, 9 is not an invariant quantity, and we can even scale it arbitrarily. 
Then, the question of the compatibility conditions is not a key issue. 

Finally, the consequences of the relation (^) lead only to two (instead of three as 
in the regular case) complex equations 

£C 6c+ l dU a b ^ Vabc = ^ (32) 

W ab CV ab - c = . (33) 

The second one is the same as the second one in the regular case. We have given their 
explicit form in the NP formalism in |Appendix A] [Eqs. (|A . 3 2|jA733| ) ] . 



4. General framework 

We have seen how the introduction of the Papapetrou field allows us to gain new insights 
into the study of the consequences of an isometry in spacetime local structure. This was 
already considered in |6|, [jj for the case of vacuum spacetimes. There, the Papapetrou 
field was used to set up a classification of these spacetimes and to study the possible 
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relations between the existence and structure of a Killing symmetry and the algebraic 
classification of spacetimes, the Petrov classification [JT3], 0] . Moreover, a new formalism 
was set up, which provides control both on the structure of the Killing symmetry and 
on the algebraic structure of the spacetime. The way this formalism works was shown 



in and a number of results were derived (see for the correction of one of the 
results). In what follows we will set up a framework which takes into account all the 
information we have presented in the previous sections, and which extends the formalism 
developed in @ for vacuum spacetimes to any spacetime, without restrictions on the 
matter content. 

The main idea is to extend the Newman- Penrose (NP) formalism by adding 
quantities, and their corresponding equations, that describe the Killing symmetry and 
its associated Papapetrou field. Briefly speaking, the Newman-Penrose formalism || is 
a tetrad formalism where all the tetrad vectors are null vectors. It has reveal itself as a 
very powerful tool in many different tasks (e.g., exact solutions, study of gravitational 
radiation, etc.), an in particular, using it we get control on the algebraic structure of 
the spacetime. The variables and equations in this formalism are (see 0, |23], ^4] for 
definitions and other details): (i) The components of a NP basis {z a b ) = (k b , £ b , m b , rh b ) 
in a coordinate system {x a }. The equations for them come from the definition of 
the associated connection. We obtain them by applying the commutators of the NP 
basis vectors (equations (7.55)-(7.58) in 0) to the coordinate system {x a }. (ii) The 
components of the connection, known as spin coefficients. They are 12 complex scalars 
named as: (k, a, p, e, u, A, fi, 7, r, ir, a, (3). The equations they obey come 
from the integrability conditions for the variables in (i), the so-called NP equations 
(equations (7.28)-(7.45) in 0), which are just the expressions of the components, in 
the NP basis, of the Riemann tensor in terms of the complex connection, (iii) Riemann 
curvature variables. These are divided into Weyl tensor components, described in the NP 
formalism by five complex scalars ^ a (A = 0, . . . , 4); components of the traceless Ricci 
tensor, given by complex scalars $jy (X,Y = 0, 1,2) satisfying $ xy = Qyx] and the 
scalar curvature, A. The equations for these quantities are the second Bianchi identities 
(equations (7.61)-(7.71) in 0, |25|), which are at the same time integrability conditions 
for the NP equations. With these equations we get a closed system of equations for 
the whole set of variables. Finally, when we get a specific description of the energy- 
momentum content through a set of matter fields, we can express the quantities $xy 
and A in terms of these fields and to add the equations for the matter fields. A good 
example is the case of electromagnetic fields (see, e.g., 0). 

Our formalism is an extension of the NP formalism based on the following two 
ideas: (A) To adapt the NP formalism as much as possible to the existence of a Killing 
symmetry. (B) To extend the NP formalism by adding new variables and equations 
characterizing the Killing symmetry. 

The best way of implementing the point (A) is through the choice of the NP basis 
{k,£,m,fh} . We will choose the NP basis to be adapted to the algebraic structure 
of the Killing symmetry. That is, in the case regular case we will choose a NP basis in 
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which the Papapetrou field F takes the canonical form (j2~ID , and in the singular case 
we choose a NP basis in which F has the canonical form (|22|). Then, in our framework 
all the quantities and equations of the NP formalism will be written in such an adapted 
NP basis. With regard to point (B), we extend the NP formalism by adding, to the 
sets described above [(i)-(iii)], the following variables and equations that describe the 
Killing symmetry: (iv) The components of the KVF in the chosen adapted NP basis. 
They are introduced through the following expressions 

Zk = k a £a, 6=^6*, £m = m°£a, U = fff^a , (34) 

where Ck and £j are real scalars, and £ m and Cm are complex and related by Cm = Cm ■ 
Then, 

£ = -£,ik - Cd + InJn + Cm-fn . (35) 

The equations for the variables (Ck,ChCm) are obtained from the definition of the 
Papapetrou field (H). Their form once projected onto an adapted NP basis is shown in 
equations (|A.1| - [A.10D of |Appendix A.l| . (v) The components of the Papapetrou field in 
the adapted NP basis. In the regular case these components are completely described 
by the complex quantity T = — (jz£ + iffi), and in the singular case by the complex scalar 
9 . The equations for these quantities come from the projection of Maxwell equations (|4]) 
onto the adapted NP basis. We have given them in [Appendix A.2| . The equations for the 
regular case are ( A.11|4ATT4]) , and the equations for the singular case are (|A.15| - [A~T8|) . 

Another set of equations that we have to take into account to complete the 
formalism are those coming from Eqs. (§|,0). As we have seen before, in the regular 
case they lead to 3 complex algebraic equations [Eqs. ( |28| - pOD 1. When we express the 
equations in the NP formalism, using a basis adapted to the Papapetrou field, they 
become equations involving the components of the Killing, the spin coefficients and 
components of the Ricci tensor [see Eqs. (|A.29| - |A~3~lD in Appendix A| . In the singular 
case we only get two complex equations [Eqs. (|32]j33| )1, and one of them, Eq. (|32f) , 
coincides with one of the regular case, namely Eq. (^). Their expressions in the NP 
formalism are similar to those in the regular case, with the only difference that now 
they involve derivatives of 9 [see Eqs. ( |A.32| , |A~33| ) in | Appendix A|| . 

So far we have presented the formalism, now let us see what is the best way of using 
it. Following the discussion in section [|, the integrability conditions for the components 
of the KVF are a good starting point since they lead to an expression for the Weyl 
tensor [see Eqs. (n,lS)]. That means that we can express the complex scalars ^> a 
(A = 0, . . . , 4) in terms of the rest of quantities of the formalism. In |Appendix A.3j we 
give these expressions for the regular [Eqs. ( |A.19| - [A~2"3"|) 1 and singular [Eqs. ( A.24JA~28| )1 
cases. 

Remarkably, the expressions for the regular case are algebraic in the spin 
coefficients, the eigenvalues of the Papapetrou field $ + i0, the KVF components, 
and the components of the Ricci tensor. The reason for that comes from the Maxwell 
equations, which are equivalent to the equations (|4]). Then, the quantity F ab . c can 
be written in terms of F a b, the connection and the Ricci tensor, hence the algebraic 
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expressions for In the singular case some derivatives of the complex scalar 9 appear 
in the Weyl tensor but, as we have mentioned before, we can arbitrarily scale this 
quantity by using the remaining freedom in the choice of the NP basis. 

The main consequence of having an expression for the Weyl tensor is that we do 
not need the equations for it. Then, from the second Bianchi identities we only need to 
consider the subset corresponding to the energy-momentum conservation equations 

G ab . b = T ab . b = 0. (36) 

Then, the situation after introducing the information coming from the integrability 
conditions for the components of the KVF is a simplified framework, where we can 
forget about the components of the Weyl tensor and their equations. 

Apart from the equations we have listed above, we must take into account the 
integrability/compatibility conditions for the components of the Papapetrou field F ab , 
and the consequences of the fact that F ab is invariant under the 1-parameter family 
of diffeomorphisms generated by the KVF £ [Eqs. (0,0)]- The first group of equations 
have been considered in section |3|. In the particular case when the Papapetrou field is 
regular and dh = 0, we have seen that they lead to five independent equations only for 
the connection, the spin coefficients, which complement the Newman-Penrose equations. 
We have given them in the |Appendix A| [Eqs. ( |A.34| - rA~3~8|) 1. In the general regular case 
the integrability conditions are (|27|). When the Papapetrou field is singular, we can only 
talk about compatibility conditions. 

Finally, in the case that we know the particular description for the energy- 
momentum content, that is, the matter fields and the field equations they satisfy, we can 
adopt the following point of view. First, we can add the matter fields as new variables of 
the formalism, and the field equations they obey as equations of the formalism. Then, 
we will get an expression, via the energy-momentum tensor, for the components of the 
Ricci tensor, $xy , A, in terms of the matter fields and, probably, of their derivatives. 
The advantage of this is that now we can ignore the quantities $xy and the equations 
that they satisfy, namely the contracted second Bianchi identities (|36"D . The last point 
would be to consider the integrability/compatibility conditions for the equations of the 
matter fields, in a similar way as we have done with the equations for the Papapetrou 
field in the previous section. In the next section, we illustrate this alternative procedure 
with the case of spacetimes with a gravitating electromagnetic field. 



5. Particular case: gravitating electromagnetic field 

To illustrate the framework we have described above we are going to consider the 
particular case of spacetimes whose energy-momentum content corresponds to an 
electromagnetic field, i.e., we have a gravitating electromagnetic field. We will call 
it H ab , to distinguish it from the Papapetrou field F ab . The self-dual field H ab will 
satisfy Maxwell equations 

H ab , b = S a , (37) 
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where S a are the electromagnetic sources. In general, H ab will not be aligned with the 
Papapetrou field, so in the NP basis adapted to the Papapetrou field it will take the 
following generic form 

H ab = %U ab + ^W ab + $ 2 Ki , (38) 

where $x (X = 0,1,2) are the components of H ab in such a basis. The energy- 
momentum tensor of the spacetime will be given by 

T ab = \H a c H cb g ab T ab = 0. (39) 

Therefore, the components of the Ricci tensor in terms of the components of H ab are 
(see, e.g., [§) 

$ xy = $ x $ y , A = . (40) 

The NP equations for $x can be found, e.g., in [Eqs. (7.46)- (7.49)]. Following the 
framework described above, we have to include the quantities $x and their equations, 



but then we can forget about which are given by (pED"D , and their equations, the 

contracted second Bianchi identities (p6|). So we have to concentrate on the Maxwell 
equations for the gravitating electromagnetic field. If we insert (|3~8"D into Maxwell 



equations (|3~7D we obtain 

$i >a = $i/i„ + q a , (41) 
q a = W ab S b + W ac ($ U bc + § 2 V hc ), b , 

where h a is given in (P5|). Again we can follow |T^] to study the integrability conditions 
for (pi), and the classification used for the case of the Papapetrou field depending on 
whether h[ ajb ] vanishes or not. Here, it is important to remark that this 1-form h a is the 
same for both the Papapetrou field F ab and the gravitating electromagnetic field H ab . 
This means we can analyze simultaneously the integrability conditions for F ab (in the 
regular case) and H ab . Then, when dh = the integrability conditions are given by 



(^q) and 

q[a,b] + Qlah] = . 

The condition dh = takes place in the case of Einstein-Maxwell spacetimes (S a = 0) 
where F ab and H ab are regular and aligned (see the particular case below in the 
subsection |5.1|) . 

In the generic case, when /i^j ^ 0, the integrability conditions are given by ([2~7j) 

and 

Qabh{ c ,d[;e ~ fya,£>] [^cd;e ~ h-c^de + <jch{d,e\] = , 

where 

^ab = Q[a,b] + Q[ a hb] ■ 
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5.1. A particular case: Einstein- Maxwell spacetimes with alignment between the 
Papapetrou and the electromagnetic field 

Let us now consider a very particular case, in which the gravitating electromagnetic 
field is source-free (S a = 0) and it is aligned with the Papapetrou field in the following 

sense 

H ab = (F ab , (42) 

where £ is an arbitrary complex scalar. In the case of a regular Papapetrou field, 
will have two principal directions which coincide with those of F a b. Moreover, we will 
have $o — ^2 — in (j3§D , which means that now we have q a = 0. Therefore, from (fiC|) 
it follows that 

$00 = $01 = $02 = $12 = $22 = = A . 
Or in other words, the only non-zero component of the Ricci tensor is $n = 
Comparing the equations for the Papapetrou and the electromagnetic fields we get an 
equation for the proportionality factor £ 

dlogC = -T V 

The integrability conditions for ( are the same as those for T, i.e. Eqs. (|26|) . 
corresponding to the case dh = 0. The interesting characteristic of this particular 
case is that if we look at the expressions for the Weyl tensor components, Eqs. ( |A.19| - 



A.23|) , we see that they are formally the same as those for the vacuum case • Also the 

NP equations coming from (||), Eqs.( |A.29|4A.31|) , have the same form as in the vacuum 

case. Finally, the particular case where the proportional factor is constant, d£ = 0, is 

empty since it implies either T = or £ = 0. 

An interesting example in which we find this kind of alignment is the case of the 

Kerr-Newman metric |26|, |2[], which in Boyer-Linquist coordinates reads 

A sin 2 6 p 2 

ds 2 = — -(dt - a sin 2 9d<p) 2 H — [(r 2 + a 2 )dp - adt} 2 + ^-dr 2 + g 2 d6 2 , 

g 2 q A 

where A = r 2 — 2Mr + a 2 + Q 2 and g 2 = r 2 + a 2 cos 2 9, being M, a and Q constants 

denoting respectively the mass, angular momentum, and charge of the Kerr-Newman 

black hole. It is easy to see that the Papapetrou field associated with the timelike KVF 

£t = d/dt has the following form 

a(2Mr - Q 2 ) sin(2#) ,„ , ^ . 2 2 , , . 
Q 

The Kerr-Newman metric is a solution of the Einstein-Maxwell system of field equations, 
and the electromagnetic field generating the energy-momentum distribution is given 
by 

H = - V2 — A (r 2 - a 2 cos 2 0) (dt - a sin 2 6dip) A dr 
Q 

-V2 sm(26)d0 A (adt - (r 2 + a 2 )dcp) , (44) 



\M(r 2 - a 2 cos 2 9) - rQ 2 ](dt - asm 2 9dip) A dr 
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which can be generated by the following vector potential 

A = -V2 r jt(dt - a sin 2 9d<p) = -^^«t + dt) . 



From the expressions for F [Eq. (|43|) 1 and H [Eq. (0)] we can construct the respective 
self-dual 2-forms, and from ( ff2"D we can find £. The result is 
Q r + ia cos 
~ V2M(r + iacos9) - Q 2 ' 

Since -F and H are aligned in the sense of (f42|) they have the same principal null 
directions. These directions, which also coincide with the two multiple directions of the 
spacetime (of the Weyl tensor), are given by the following vector fields 

A asin 2 6A J _ A Jjt , asin 2 #A J 
K = — -dt + dr-\ dcp , L = — -dt - dr H dip . 

Q z Q A Q z Q z 

Note that these vector fields are not normalized, i.e. K a L a ^ —1. From them we can 
compute the complex scalar T containing the eigenvalues ^ and j8 [see Eq. (|2T|)1: 

(r + ia cos 9) 



T = -2- 



M(r + ia cos 9) - Q z 



Using all this information we can finally get q and h\ 

Q 2 

q = —2—d(r + ia cos 0) , h = —2d \n(r — ia cos 6) . 

Q 

As expected, we have dh = 0. 

In the case of a singular Papapetrou field, H a b will be also singular and with the 
same principal direction. In this case we have $ = = in (p8|), hence we find that 

$oo = $01 = $02 = $n = $12 = = A . 

Then, the only non-zero component of the Ricci tensor is $22 = $2$2- From the 
expressions for the {A = 0, ...,4) given in the |Appendix Aj we have 

^0 = #i = ^2 = , 

which means that the spacetime is not only algebraically special, as the Mariot-Robinson 
theorem [27, |28|] tells us, but also that it must be either of the Petrov type III or N. 
In other words, the only principal null direction, k, must have at least multiplicity 
three. This generalizes the result found in for vacuum spacetimes. Here, when the 
proportional factor ( is constant, one can see that this leads to the additional condition 
p — 0, which means that we are within the Kundt class of spacetimes [p9fl. 



6. Remarks and discussion 



In this paper we have discussed the consequences of a symmetry on the local structure 
of the spacetime. We have adopted the point of view introduced in || 0], where the 
Papapetrou field plays a central role. We have also generalized the formalism introduced 
in for vacuum spacetimes to the general case. In our discussion we have shown that 
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the advantages of this type of formalism lie on the consideration of the following two 
essential points: first, that the spacetimes of General Relativity are four-dimensional, 
which we have shown leads to an expression for the components of the Weyl tensor 
through the integrability conditions for the components of the KVF. In this line, the 
formalism we have presented is based on the Newman-Penrose formalism, in which the 
fact that the spacetime is four-dimensional is present. The second important point is 
that we have written the NP equations in a basis adapted to the algebraic structure of the 
Papapetrou field, and we have added variables, with their corresponding equations, that 
describe the KVF as well as its associated Papapetrou field. By one hand this simplifies 
the equations, specially if we want to study cases in which there are alignments between 
the principal directions of the Weyl tensor and the Papapetrou field (see [^|, |14|)- O n the 
other hand, these new variables provide a new point of view to study the integrability 
conditions. Finally, using the Papapetrou field we introduce information about the 
particular characteristics of the symmetry, for instance, when we decide whether a KVF 
has a regular or a singular Papapetrou field. 

We have illustrated the formalism in the case of spacetimes with a gravitating 
electromagnetic field, where we have described how the formalism should be applied 
and we have produced some new results. The fact that this formalism can be applied 
to any spacetime with a Killing symmetry means that we can use it with a number 
of interesting physical situations, and in this sense it could be used to complement or 
in combination with other methods of study, like perturbative schemes or numerical 
computations of spacetimes. 

Finally, it is important to remark that here we have dealt only with the case of 
Killing symmetries however, the structure of the formalism and the way in which it 
has been set up, can be easily translated to other kinds of symmetries. Actually, there 
are some other works in which this has already been done. The case of homotheties 
in vacuum spacetimes has been studied in ||15| , and in |16|] it was applied to conformal 
Killing symmetries. 
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Appendix A. Formalism for spacetimes with a Killing symmetry. NP 
equations 

In this Appendix we give the equations of the general framework present above using 
the Newman-Penrose formalism. 
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Appendix A.l. Equations for the components of the KVF 

The equations for the components of the KVF [see Eqs. (|34|,|35])] are obtained by 
projecting Eq. (f|) onto an adapted NP basis. The result is the following (they are 
the same as those found in @ for the vacuum case): 



D£k - (e + e)6fc + K^m + «£m = 


0, 


(A.l) 


A£ fc - (7 + 7)6; + r£ m + r£ m = 


= if*> 


(A.2) 


- (« + P)£k + Pirn + Crlm = 


0, 


(A.3) 


+ (e + e)& - ir$ m - 7r£ m = 




(A.4) 


+ (7 + 7)6 - Km ~ Km = 


0, 


(A.5) 


5& + (a + (3)£i - - \£ m = 




(A.6) 


Dim - (e - e)^m - Kk + = 


0, 


(A.7) 


A£ TO - (7 - l)im ~ Kk + Til = 




(A.8) 


Sim + {a- f3)im ~ Mk + o~ii = 


0, 


(A.9) 


sim - (a- (3)im - fiik + pii = 




(A.10) 



where (D, A, 5, 6) are the directional derivatives along the NP basis vectors 

D = k a d a , A = e a d a , 5 = m a d a , 5 = m a d a . 

These equations include at the same time both the regular and singular cases. The 
regular case is obtained by taking 9 = and for the singular case we need to take 
ft = 18 = 0. 

Appendix A.2. Equations for the components of the Papapetrou field 

They are the projection of the Maxwell equations (|) onto the adapted NP basis. 
In the regular case they are equations for the complex quantities $ + actually 



equations ( P^ ) : 

D(jt + = 2p({& + ifS) - 4 {&($n - 3A) + 6$oo - £ m <l> 01 - £ m $ 01 } , (A.ll) 

A(*fc + i0) = - 2^(f& + 1 (3)+ 4 {i k $ 22 + i^n - 3A) - im$i2 - im^n} , (A.12) 

6{fl + i0) =2r(^ + 2^)-4{a$i 2 + 6$oi-em( < fii + 3A)-e m < fo 2 } , (A.13) 

5(fi + ip) = -27r(^ + i^)+4{c fe $i2+^oi-^02-L($n+3A)} . (A.14) 



Then, we have an expression for all the directional derivatives of (Ai and In the 

singular case, Maxwell equations reduce to equations fl3~l|) for 9, which once project 
onto the adapted NP basis are 

D9 ={p- 2e)9 + 2 {&$ 12 + - £ m $02 - £m($n + 3A)} , (A.15) 

59 = (r - 2/3)0 + 2 {£ fc $ 2 2 + £,($n - 3A) - £ w $ 12 - £^12} , (A.16) 

\k9 = &($ii - 3A) + e^oo - im%i - L%i , (A.17) 

\o9 = i k $ 12 + - &„($n + 3A) - £ m $ 02 • (A.18) 
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In this case, we only get two directional derivatives of 9, in the directions of k and 
m. The other two equations are algebraic relations, from where we can recover the 



Goldberg-Sachs and Mariot-Robinson theorems [[3(], 27, 28] in the particular case of 
spacetimes with a KVF. 

Appendix A. 3. Integrability conditions for the components of the KVF 

The integrability conditions for the components of £ lead to an expression for the Weyl 
tensor [Eqs. ( |TT| , |T~8|)] . From it we can compute the Weyl complex scalars ^a ■ In the 
regular case, we get the following expressions: 

*o = ^^{<m - atk) - §(d$oo - + 02) , (A.19) 

*i = - *U) - ^ {aCm^oo - (66 + U£ m )®oi + 66^02} , (A.20) 

* 2 = - tL) - I {tf%o - 2&£ m $ i + - 2A , (A.21) 

*s = ^^-Hk - Km) - § {teA - (66 + 6»£m)$i2 + 66^22} , (A.22) 

* 4 = ^(4 - A6) - |te 2 $02 - 266^12 + e 2 $ 22 ) • (A.23) 

Here we have used that N = — 2£fc6 + 2£ m £ m . In the singular case, the complex scalars 
*f?A can be written as follows: 

*o = ~^&oo - 2£ fe £ n $ i + £^02) , (A.24) 

1 = -— {66^00 - (66 + 6*60*01 + 66n$0 2 } , (A. 25) 

^2 = ~ {tf%o - 266^oi + &02} - 2A , (A.26) 
^3 = ~{f>£i ~ rU)6 - ^(6 2< Soi - Hilm*ix + &12) , (A.27) 
M> 4 = 1 {e m (A + 2 7 ) - 6(5 + 2a)} - ^(6 2 $02 - 266^12 + 6^22) ■ (A.28) 

Appendix A. 4- Invariance of F a b 

From the Ricci identities we get the equations for the Papapetrou field, the components 
of the Weyl tensor and some extra relations, equations ( p8[ - |30l) in the regular case and 
equations (|32| , |33f) in the singular case. These relations come from the invariance of F a & 
[Eq. (H])]. In the NP formalism, the equations for the regular case are: 

+ <i ~ Km -(i£ m = 0, (A.29) 
t~6 - ptm + «6 - o-£ m = , (A.30) 

A*6 - /°6 - ^61 + T£ m = 

-rt-l,&22 - 6^00 - 26Cm$12 + 26? m $01 + d$20 " ^$0 2 ) • (A.31) 
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And for the singular case: 

r£k ~ Pim + k£i ~ cr^ m = , (A.32) 
&(A + 2 7 )# - £ m (5 + 2a)9 + 06 - rU)6 = 

- 2& 2 $oi + 2&U%2 ~ 2(&£, + CmL)$i2 + - &12 + 26^ m $ 22 • (A.33) 

Appendix A. 5. NP form of the equations dh = 

Within the regular particular situation of interest corresponds to the subcase 

defined by the condition dh = 0. These equations can written in the NP formalism in 
the following form: 

(D + e + e)fji + (A - 7 - f)p - nn + rf = , (A. 34) 

(5 + 7r - a - (3)p - (D - p - e + e)r + fin - na = , (A.35) 

(5 - a - (3)p + (D + e - e)vr + fiR + to = , (A.36) 

{5 + a + p)n + (A - 7 + 7)7 - pv - nX = , (A.37) 

(A + p, + 7 - 7)tt - (5 - f + a + (3)p - t\ + pv = . (A.38) 
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